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1. Introduction 

Non-perturbative investigations of Higgs-Yukawa models using the lattice regularization be- 
came subject of many investigations in the early 1990's [fil 0, 51 ||, ||, §]. These lattice studies were 
initially motivated by the interest in a better understanding of the fermion mass generation via the 
Higgs mechanism on a non-perturbative level and, in particular, in the determination of bounds on 
the Yukawa couplings translating into bounds on the Higgs boson mass and the - at that time not yet 
discovered - top quark mass. Furthermore, the fixed point structure of the theory received special 
attention due to the question whether besides the Gaussian also a non-trivial fixed point might ex- 
ist. However, these investigations were blocked, since the influence of unwanted fermion doublers 
could not successfully be suppressed. Moreover, the models of these studies suffered the lack of 
chiral symmetry. The latter, however, would be indispensable for a consistent lattice regularization 
of chiral gauge theories like, for example, the standard model of electroweak interactions. 

Here, we follow the proposition of Liischer ^ for a chiral invariant lattice Higgs-Yukawa 
model based on the Neuberger overlap operator [Q]. As a first step we begin with an analytical in- 
vestigation of its phase structure by means of 1 /./Vy-expansions following [^, [T^] and the references 
therein. We derive an expression for the effective potential at tree-level and present our preliminary 
results for the corresponding phase diagram. 

In chapter || we briefly describe the considered model, before we present the large Nf results 
at tree-level for small Yukawa couplings y N <x 1/ y/Nf in chapter [3| In the following chapter |] we 
discuss a different large A^-limit which becomes valid for non-vanishing Yukawa couplings. We 
then end with a short outlook. 

2. The model 

The chiral invariant Higgs-Yukawa model, which we consider here, contains one four-component, 
real Higgs field <I> and Nf fermion doublets represented by eight-component spinors \\r- 1 ' , xjr- 1 ) with 
i = l,...,Nf. Furthermore, there are also Nf auxiliary fermionic doublets %^ only introduced 
to construct a chiral invariant Yukawa interaction term. The partition function can then be written 

as 

exv^-Sv-Sf-Sy) (2.1) 

where the total action is decomposed into the Higgs action 5$, the kinetic fermion action Sp n , and 
the Yukawa coupling term Sy. It should be stressed that no gauge fields are considered here. 

The kinetic fermion action describes the propagation of the physical fermion fields y/w,^ 1 ) 
in the usual way according to 

<,kin _ V V vr® ® {ov] w W _9ny (i) 1 v W 
— 7^ 7. yn -^n,m ym ^-HAn ^njnfon 

i=\n,m 

where the coordinates n,mas well as all field variables and coupling constants are given in lattice 
units throughout this paper. The (doublet) Dirac operator f^( ov ) = f^( ov ) (g) Qj( ov ) i s given by the 



Z= I D<$>Y[ D V {i) D ¥ (i) D X [i) D X {i) 
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Neuberger overlap operator Q)^ ov \ which is related to the Wilson operator = y^V^ — rV^V^ 
by 



A } 



^) =pil + -—K A = ^ w >-p, p>\ (2.2) 

where denotes the symmetric difference quotient. It is well known that the eigenvalues v ± (p) 
of £F( ov ) with Im[v ± (/?)] ^ form a circle in the complex plane, the radius of which is given by the 
parameter p. In momentum space with p £ = [—71, 7i]® 4 these eigenvalues are explicitly given 
by 

v ± {p )= p + p . ^^+ r f -P p,=sm(p,), ^=2sin(^). (2.3) 
sf p 1 + (rp l - p) 1 v £ ' 

The auxiliary fields X on the other hand do not propagate. They do not have a direct physi- 
cal interpretation at all and their contribution to the action is only introduced to establish chiral 
symmetry. 

The Higgs field couples to the fermions according to the Yukawa coupling term 



N f 

(0 , ~(<'h 



n.mi=\ 



(1-75) (1 + 75) +' 



(¥m+Xm) (2-4) 



where y N denotes the Yukawa coupling constant and B n m will be referred to as Yukawa coupling 
matrix. Here the Higgs field <!>„ is rewritten as a quaternionic, 2x2 matrix (j> n = <&„1 — i<t> J n Tj 
(Zj: Pauli matrices) acting on the flavor index. Due to the chiral character of this model, left- and 
right-handed fermions couple differently to the Higgs field, as can be seen from the appearance of 
the projectors (l±)§)/2in the Yukawa term. 

Finally, we use a slightly unusual notation for the Higgs action given by 

s* = -k*J>2 [$ n+A +$„_ A ] -Nff (2.5) 

n,n n n 

where Kjv denotes the hopping parameter and Xn is the quartic coupling. The usual notation is 
reobtained by a trivial rescaling of the coupling constants and the Higgs field. The model then 
obeys an exact, but lattice modified chiral symmetry recovering the actual chiral symmetry in the 
continuum limit [0]. 

For the further analytical treatment of this model the fermionic degrees of freedom are in- 
tegrated out leading to the effective action S e ff[&], which can be written in terms of fermionic 
determinants yielding 

S eff [<t>] = S*[$\ -N f ■ log [det (y N B&W - 2p¥ ov) - 2py N B* 



(2.6) 



3. Large A//-limit for small 



We now consider the limit of infinite fermion number Nf — > °° with the coupling constants 
scaling according to 

yN = — , 5>n = const Xn = #, = const % = kfr, = const (3.1) 



N f > 
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which justifies a semi-classical approach to the investigated model. We will therefore directly 
evaluate the effective action ( |2.6| ), which is at least possible for the constant and the staggered 
mode of the Higgs field. Since the Higgs field scales proportional to y/Nf in the large Nf -limit, we 
apply the ansatz 

$n = $- y/Nf- (m+s-i-ifi*} (3.2) 

where $£R 4 with |4>| = 1 denotes a constant 4-dimensional unit vector. We will refer to m, s as 
magnetization and staggered magnetization, respectively. For the actual evaluation of the effective 
action we rewrite (2.6), neglecting all constant terms independent of <J> in the following, as 

-l 



S eff [<i>} = S*[<!>]-N f -log 



det 1 



2p 



■B 



(3.3) 



For the effective potential V e ff(m,s) at tree-level one then finally finds 



—V eff (m,s) 



in 



d 4 p 
(2kT 



log 



+ m 2 + s 2 + X h 



4 4 r 2 2 
m + S + 0/77 s 



(m 2 + s 2 ) 



2pJ 



(m 2 -s 2 ) 



\v(p)-2p\ \v{p)-2p\ 



|V(P)| 



\y(P)\ 



2p 



v(p)-2p\ \v(p)-2p\ 



Mp)\ 



\v(P)\ 



(3.4) 



with the abbreviations pn = + %. 

The phase diagram can then be explored numerically by searching for the absolute minima of 
the effective potential V e ff(m,s) with respect to m and s. In general, four types of solutions can 
be obtained. These are a symmetric (SYM: m = 0, s = 0), a ferromagnetic (FM: m / 0, s = 0), 
an antiferromagnetic (AFM: m = 0, s / 0), and a ferrimagnetic phase (FI: m ^ 0, s ^ 0). The 
corresponding phase diagram for the choice = 0.1 is shown in Figure [I]. Here, we only present 
a qualitative plot of the phase structure, since the obtained results are preliminary. However, the 
diagram reveals a rich phase structure. Especially, there is a ferrimagnetic phase, which was also 
observed in Monte Carlo studies of earlier, not chiral invariant Higgs-Yukawa models, for example 
in M. 



4. Large iV/-limit for finite 



We now consider a different limit of large fermion numbers where the coupling constants scale 
according to 



yN = yN,5>N = const X N = , X N = const k n = ^-, k N = const. 



(4.1) 



Again, the Higgs field scales proportional to WNf allowing to expand the effective action in powers 
of 1/ (yjvl^l). In this large Nf -limit the power series expansion 



S eff [0] = —Nf f Tr log [B] - £ - (JL\ Tr ^ - 2p) B 



(4.2) 
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K N 




Figure 1: Qualitative phase diagram with respect to the Yukawa coupling y# and the hopping parameter Kjy 
for a constant quartic coupling Xn = 0.1. 



of the logarithm in equation (3.3) can therefore be cut off after the first non-vanishing term. The 
first summand (k = 1) is identical to zero and the second term (k = 2) is the first non- vanishing 
contribution. Cutting off the power series at k = 2, the model then becomes an effective, 0(4)- 
symmetric spin model 



/ (?n) 2 <t> M <t> v 



yN n,i 



<T> 12 



(4.3) 



where Ax = n — m and the non-local coupling matrix K^ v (Ax) = 5^ ]V • |r(Ax)| with |.| denoting 
the 4-vector norm is explicitly given by the momentum integral 

v+(p) 



RbT^(Ax) 



d 4 p 



jpAx 



pel 



(271) 



v+(p)-2p ^2' 



(4.4) 



which can be computed numerically. One then finds that the square-norm of the coupling matrix, 
|£(Ax)| = |r(Ax)| 2 , decays exponentially with increasing distance \ Ax\ as shown in Figure ^[ In a 
field-theoretical sense the effective spin model therefore remains a locally interacting model. 

For the evaluation of the corresponding phase diagram of this effective spin model we make 
the ansatz 

&n = y/Nf ■ <P0 ' c« 5 Cn G 1R 4 , \a„\ = 1, R 3 (po = const (4.5) 

for the Higgs field. Considering only the leading power in l/Nf of the tree-level effective action 
S e ff[0], the amplitude <po can be fixed according to 







-A-\ + \+2l N -((pl-\) 



(4.6) 



Including the next to leading order terms in 1 /Nf, the model effectively becomes an 0(4)-symmetric 
non-linear a -model 



S e f f [*]=-£ K% ■ °n ■ <5,n , With 



run 
,2 



I6p z ±4 
<{m = ^2 ' l r (Ax) I 2 + % • 9o ■ L ^ m+ft , 



(4.7) 
(4.8) 



n=±i 
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the phase structure of which is again accessable to an 1 /./V-expansion. Here, N denotes the number 
of components of the Higgs field <1>, i.e. eventually N = 4. As usual, the constraint \a n \ = 1 is 
removed by the introduction of an auxiliary field X n according to 



DAfJ[Da'] exp[-5[a,A]], 



i=\ 



S{a,X] = .££-<//. ^.a; + £A„ 



tN 



N 

1=1 



1 



(4.9) 
(4.10) 



n,mi=\ n 

Integrating out the N —\ components o 2 ,...,o N one obtains the reduced action S[c x , X] depending 
only on the fields X n and a,\ given by 

>2 



S[a l ,X] 



nm 'N n l N L 



1 



+ -(Af-l)Tr„, m log[-<# + AA, 



(4-11) 



The newly introduced parameter is necessary to perform the large Af-limit, which will be done 
for iff = tjf'N = const. Here, we are actually interested in the case tjf = 1 and N = 4. We therefore 
choose ?w = 4. Again, we consider the magnetization and a possible staggered magnetization for 
the o,] field and assume the auxiliary X n field to be constant leading to the ansatz 



o n \= m + s- (—iy andA„ = A. 
Minimizing S [a 1 , X] with respect to m, s, and X will then yield the three gap equations 








m 2 + s 2 



m ■ 



s ■ 



X 

X - 



, 2 , !6p 2 , 



-(1--) 
4 ^ 



k£&> 



y N % 

d 4 k 



-Kjv<p 2 £cos(^) - ^-^q{k) + ^ 

it n% 2 



(4.12) 

(4-13) 
(4.14) 

(4.15) 



where q(k) denotes the eigenvalues of |r(Ax)| corresponding to plane waves and is given by 



q(k) 



pet 



d 4 P v + (p) 

(27iyv+(p)-2p v+(p)-2p 



P-P 



p = k-p. 



(4.16) 



For the ferromagnetic phase, i.e. m/0 and s = 0, one obtains a self-consistent determination 
equation for m 2 > according to 

-l 

(4.17) 



0<^ = 1-|(1-I) 



d 4 k 



(2n) 



8p 2 



^<Po£ (1 -cos(* M )) + ^3 



An analogous relation holds for the antiferromagnetic phase (AFM). These relations can be treated 
numerically. The resulting phase diagram is shown in Figure ||. Again, we present a qualitative 
diagram only, since our results are preliminary. The appearance of three different phases, namely a 
symmetric, a ferromagnetic and an antiferromagnetic phase, can be observed. The presented phase 
diagram is in qualitative agreement with earlier results, using non-chiral fermions, see e.g. 
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Figure 2: Left: Square-norm of the coupling matrix K(Ax) versus the square-norm of the distance Ax in 
lattice units. Right: Qualitative phase diagram with respect to the Yukawa coupling ^> 1 and the hopping 
parameter Kjv for a constant quartic coupling \ N = 0.1. 

5. Summary and outlook 

In this paper we have studied the phase structure of a chiral invariant lattice Higgs-Yukawa 
model, originally proposed by Luscher, by means of analytic 1/A^-expansions at tree-level. This 
was possible for small and for large values of the Yukawa coupling. Symmetric, ferromagnetic, 
and antiferromagnetic phases have been observed in both regimes of the Yukawa coupling con- 
stant. Additionally, a ferrimagnetic phase was found in the first regime. A next step would be the 
comparison of our analytical results with corresponding Monte Carlo data. 
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